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Abstract 



COLOUR INK MODEL PROCESSES FOR PRINTERS 



A method of predicting colours resulting from using one ink on an ink jet printer is 
disclosed. The method comprises firstly producing a test page for the printer, comprising 
sample ink patches, each patch comprising ink dots being of a predetermined quantity of 
ink. Then, measuring each sample ink patches to produce colour space coordinates for 
each patch. Therafter values representing a total quantity of ink for each patch and colour 
space coordinates are fitted to predetermined functions using a minimisation process, 
wherein the functions are of the fonm P.,. (I -fx(a)) + f/a)\ Py (I - fy(a)) + ly fy(a)\ 
Pz 0 'fz(^)) h fz(a)\ where a represents the quantity of fx(a) , fy(a) Hndf~(a) are 
functions of a, and Pv, Py and P- are colour space components representing a paper colour, 
and /r, ly and /. are colour space components representing the ink colour. 
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COLOUR INK MODEL PROCESSES FOR PRINTERS 
Field of Invention 

The present invention relates to colour ink model processing including a method 
and apparatus for printer characterisation and, in particular, for characterising colour 
5 printers, and inversion of the colour ink model. The invention also relates to a computer 
program product including a computer readable medium having recorded thereon a 
computer program for printer such processing. 

Background 

Colour printer characterisation is the process of determining the colour a printer 

10 will produce when a certain amount of ink of the various available colours is requested by 
a given printer. In order to carry out colour printer characterisation, a relationship 
between signals input to the printer and colorimetric values for printed colours must be 
established. The relationship is generally expressed as a printer characterisation function. 
One known method of determining the printer characterisation function for a given printer 

15 is by firstly, producing a test page with a predetermined number of sample colour patches, 
secondly, measuring the colour of each colour patch and finally, interpolating among the 
measurements. However, this method is time consuming and is not very accurate. 

Another known method of determining the printer characterisation function is to 
determine a printer model, which approximates the function. The main advantage of 

20 modelling is that the printer characterisation can be carried out with a comparatively 
small number of measurements resulting in a large time saving. In determining the 
printer model, some assumptions are necessary in order to simplify the mathematics 
involved. The accuracy of the printer model will be limited by these assumptions. 

One known printer model is the Neugebauer mixing model. The Neugebauer 

25 model is used in modelling and characterising colour printers. The Neugebauer model is 
used to predict the colour of a print on a given printer, as a weighted average of the XYZ 
values of the solid overprints of the three primaries (ie. cyan (C), magenta (M) and yellow 
(Y)). The weights of each colour are determined by the relative dot area coverages of C, 
M and Y constituting the print. The dot area coverages for the digital input values are 

3D determined using a combination of direct measurement and calculation. 

As discussed above, the Neugebauer model provides the characterisation 
function from device values (C,M,Y) to colorimetric values (XYZ). However, for printer 
characterisation the inverse mapping (ie., from colorimetric values to device values) is 
required. 
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Numerical methods are used to invert a printer characterisation function, which 
is non-linear. However, a problem occurs when a printer characterisation function has 
more than three inks, as a number of different ink combinations can result in the same 
colour. 

Several methods have been proposed for optimising the Neugebauer model. One 
such method is discussed in an article by Balasubramanian et ah entitled "Optimisation of 
the spectral Neugebauer model for printer characterisation". Journal of Electronic 
Imaging, April 1999, Vol. 8(2). The method uses weighted spectral regression for 
optimising the Neugebauer primaries in order to characterise Xerographic printers. 
However, the method proposed by Balasubramanian et al does not work on ink-jet 
printers where inaccuracies are found in the optimised Neugebauer model proposed in the 
article. 

It is an object of the present invention to ameliorate one or more of the 
limitations of the methods described above. 

Summary of the Invention 

In accordance with one aspect of the present disclosure, there is provided a 
method of predicting colours resulting from using one ink on a printer, said method 
comprising the steps of: 

producing a test page for said printer, said test page comprising a plurality of 
sample ink patches, each said patch comprising a plurality of ink dots according to a 
predetermined ratio for each patch, said predetermined ratios being based on a number of 
dots printed for a corresponding patch compared to a maximum number of dots; . 

measuring each of said sample ink patches to produce a plurality of colour space 
coordinates for each said patch; and 

fitting values representing a ratio of dots for each said patch, and said colour 
space coordinates, to a plurality of predetermined functions using a minimisation process 
in order to predict said colours, wherein said predetermined functions are of the form 

Px (1 -m) + h fM \ 

Py(l'fy(a)) ^Iyfy(a)\ 

Pz 0 'fz(a)) + h f.(ah 
where a represents a ratio of dots for a particular \^^X^yfx(o),fy(a) and fz(a) are functions 
of a and are derived from a section of the curve y = //x, and wherein P^, Py and Pz are 
colour space components representing a paper colour, and I^^ ly and are colour space 
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components representing the ink colour, wherein said fiinctions fj(a)^fy(a) and fi(a) are of 
the form: 

and wherein a is a variable and k is a constant, which can be different in each oifx(a)^ 
fy(a) Bndf^(a), 

In accordance with another aspect of the present disclosure, there is provided an 
apparatus for predicting colours resulting from using one ink on a printer, said apparatus 
comprising: 

means for producing a test page for said printer, said test page comprising a 
plurality of sample ink patches, each said patch comprising a plurality of ink dots 
according to a predetermined ratio for each patch, said predetermined ratios being based 
on a number of dots printed for a corresponding patch compared to a maximum number 
of dots; 

means for measuring each of said sample ink patches to produce a plurality of 
colour space coordinates for each said patch; and 

means for fitting values representing a number of dots for each said patch and 
said colour space coordinates to a plurality of predetermined functions using a 
minimisation process in order to predict said colours, wherein said predetermined 
functions are of the form 

(1 -Ua)) + f,(ay, 
Py a 'fy(a)) + ly fy(a)\ 

P2 0 -ua)) + h m\ 

where a represents a ratio of dots for a particular T^^Xohy fx(a) , fy(a) and f2(a) are 
functions of <a and are derived from a section of the curve = //x, and wherein P^, Py and 
are colour space components representing a paper colour, and 7^, ly and Ij are colour 
space components representing the ink colour wherein said functions fx((^)^fy(^) and fzia) 
are of the form: 

and wherein a is a variable and /: is a constant, which can be different in each oifx(a)^ 
fy(a) and //a;. 
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In accordance with still another aspect of the present disclosure, there is provided 
a program for predicting colours resulting from using one ink on a printer, said program 
comprising: 

code for measuring each of a plurality of sample ink patches of a test page for 
5 said printer to produce a plurality of colour space coordinates for each said patch, each 
said patch comprising a plurality of ink dots according to a predetermined ratio for each 
patch, said predetermined ratios being based on a number of dots printed for a 
corresponding patch compared to a maximum number of dots; and 

code for fitting values representing a number of dots for each said patch and said 
10 colour space coordinates to a plurality of predetermined functions using a minimisation 
process in order to predict said colours, wherein said predetermined functions are of the 
form 

Px (1 'f.(a)) + I, Ma) \ 
Py(l'fy(a)) + Iyfy(a)\ 
15 (1 -f^ia)) + /, f,(a)\ 

where a represents a ratio of dots for a particular patch, fx(o)^ fy(a) and /^(a) are 
functions of a and are derived from a section of the curve y = 7/x, and wherein Py and 
P^ are colour space components representing a paper colour, and ly and /z are colour 
space components representing the ink colour, wherein said functions fx(o)^fy(cO and f^(a) 
20 are of the form: 

and wherein a is a variable and ft is a constant, which can be different in each of fx(a)^ 
fy(a) and f,(a). 



25 
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Brief Description of the Drawings 
A number of preferred embodiments of the present invention will now be 
described with reference to the drawings, in which: 
5 Fig. 1 is a flowchart showing a method of characterising a colour printer; 

Fig. 2 shows a cube representing the Neugebauer mixing model for a three ink 

printer; 

Fig. 3 A shows the area of one dot printed on paper; 
Fig. 3B shows the area of a subsequent dot printed on paper; 
10 Fig. 4 shows the function Y=l/X; 

Fig. 5 shows a dot with a single halo printed on paper; 

Fig. 6 shows the method of solving the multi-ink model in accordance with the 
preferred embodiment of the present invention; and 

Fig. 7 is a schematic block diagram of a general purpose computer upon which 
15 the process described herein can be practiced. 

Detailed Description including Best Mode 
Where reference is made in any one or more of the accompanying drawings to 
steps and/or features, which have the same reference numerals, those steps and/or features 
have for the purposes of this description the same function(s) or operation(s), unless the 
20 contrary intention appears. 

Fig. 1 is a flowchart showing a method 100 of characterising a colour printer 
having particular application to ink-jet printers. However, the method can be used to 
characterise any printing device that perfomfis colour gradation by placing a plurality of 
dots of the same size on a page. The process begins at step 101, where a test page (or 
25 Measurement Set) is produced for a given printer to be characterised. The test page 
includes a number of sample colour patches. To ensure accuracy, a sufficiently large 
number of sample colour patches covering the fiill colour gamut of the printer, is 
required. The ink amounts for each colorant (eg. CMYK) contributing to each sample 
colour patch is also required to be known. At the next step 103, a colorimeter is used to 
30 take measurements of the sample colour patches, in a colour space which is linear in 
power of light. The CIE (1931) XYZ colour space is preferably used. However, any 
linear colour space can be utilised. The result of the colorimeter measurements is a set of 
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XYZ values for each sample colour patch on the test page. The process continues at step 
105, where the ink amounts for each colorant of each colour sample patch and the 
corresponding measured XYZ values are fitted to a function, using a minimisation 
process such as the least squares process. 

The printer characterisation function thus formed is a multi-linear interpolation 
of an extended set of Neugebauer primaries with an adjustment. The function is derived 
from the assumption that each dot printed on a page by a printer has a small halo of 
slightly different colour. The function will be explained in more detail later in this 
document. The computation must be repeated independently for X, Y and Z. The result 
of each independent computation is a set of values for each of the extended set of 
Neugebauer primaries and values for two constants and k'i for each ink. The 
computation will be explained in more detail later in this document. 

The method 100 of characterising a printer will work on a printer with any 
number of inks. For example, the preferred method will give device values for a CMY 
printer, a CMYK printer or a six-ink (i.e., C^CTvTM YK) printer. 

The printer characterisation function of the method 100 may be conveniently 
appreciated by letting S represent the set of all inks of a printer. The usual definition of 
the set of Neugebauer primaries is the set of colours of all combinations of the elements 
of the set S. The combination which is the empty set is also included and represents the 
paper colour. The method 100 uses an augmented set which the inventor has referred to 
as a "Double Neugebauer set" and is based on a set S which contains not only the inks of 
the printer, but also the halos of the inks of the printer. The basis set S for the Double 
Neugebauer set is that such a set contains twice the number of elements (as found in the 
basis set for a standard Neugebauer set) and therefore the number of elements in the 
Double Neugebauer Set is the square of the number in the usual Neugebauer set. 

Inkjet printers must use a halftoning method to determine which dots to turn on: 
In the method 100, the proportion of dots which are turned on is represented as (ay,...,a„) 
where n is the number of inks on the printer. 

The X component of the XYZ colour of a patch printed on an ink jet printer with 
n inks which have been printed with levels {a/,...,a„) for each of the inks is given by 
inserting the computed values for X for each of the Double Neugebauer primaries, and the 
values for the constants ki and into the following printer characterisation function (ie: 
"Multi-Ink Model): 
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<^ = S C Jn/.e^? /.(a,): 1 -/(«,-)] [n ^. ^ : 1 - (I) 

p^K L '=0 J L '=0 J 

where: 

K is the Double Neugebauer set of colours 

5 Cp is the X, Y or Z component of the XYZ colour of the primary p; 

/, € p is true if the i**" ink is included in the Neugebauer primary p; 

Hi e p is true if the i*** ink's halo is included in the Neugebauer primary p; and 

^^^^ = -TT^-rT' (2) 

gifaj = ki'ar, and 
10 the expression q ? a : b takes the value a if ^ is true and b 'xfq 'xs false. 

Similarly, the Y and Z components can be computed and will give rise to 
different values for A:, and ki\ 

The method 100 of characterising a printer will now be explained in further 
detail including the manner in which the Double Neugebauer Set is derived. 
15 When an ink is printed on a page as dots which have consistent colour, the 

resultant perceived colour can be expressed as (X,Y,Z), where; 

X = jcIX(>.)I(X)dX + {\-x)\X{X)V{X)dX, 
Y^x\ Y{X)l{X)AX + (1 -;c; J Y(X)P(X)dX, 

Z = x\Z{X)\{X)dX + {\-x)]Z{X)?{X)dX, (3) 
20 and where, x is the proportion of paper that is covered with ink; 

X(X), Y{X) and Z{X) are the CIE 1931 Standard Observer sensitivity 
functions; 

I(X) is the ink reflectance spectrum for the given illuminant; 
P(X) is the paper's reflectance spectrum for the given illuminant; and 
25 X, Y and Z are the resultants components of the colour in the CIE 1931 

XYZ colour representation. 
Notably, J X(X)P(X,)dX, J Y{X)?{X)dX , and \ Z{X)?{X)AX are all measurable as the 
colour of the paper for the given illuminant. This can be expressed as: 

P ^(Xp.Yp.Zp) (4) 
30 Similarly to the single ink case discussed above, for a CMYK printer the CIE 

XYZ values can be measured for each colourant and expressed as follows: 

C =(Xc. Yc^Zc); 
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Y = Yy, Z^O; and 

K = (Xj, Yi, z,y 

Qt M_y Y. and represent the colours of each ink on a given paper type for a 
5 given illuminant. 

Similarly, the colours of ink combinations can be defined. For example, let CM 
represent the colour of the paper covered by both Cyan and Magenta ink. 

For two inks (eg., cyan and magenta), the colour of the two inks printed as two 
layers on a given paper type, for a given illuminant, can be expressed as follows: 
10 0 'Cj (I -m)P -^(1 -c)mM+c(} - m) C-^cm CM (5) 

where: c is the proportion of the page covered by cyan ink; and 
m is the proportion of the page covered by magenta ink. 

The above expression (5) assumes that the ink dots are randomly positioned over 
the page. For more inks the expression (5) can be extended. 
15 Fig. 2 shows a cube 200 representing the Neugebauer primaries for three inks 

cyan, magenta and yellow. Expression (5) can be extended to correspond to a tri-linear 
interpolation in the cube 200, as shown in Fig. 2, where the colours measured are in the 
comers. 

To predict the colour of a patch of a single ink a non-linear function of the 
20 proportion of dots printed must be used. For an ink-jet printer the inventor has observed 
that the placing of ink on a given paper has some specific characteristics. Firstly, because 
the dots of ink are all of the same size, dot gain is largely irrelevant to the model once the 
effective dot size is determined. Secondly, the dots are placed in a grid pattern. Fig. 3A 
shows a representation of the area covered by one dot 301 printed on paper. The dots in 
25 broken lines (eg. 304, 305) represent the area covered by subsequent dots printed on the 
paper. Therefore, the first dot printed covers a full dot area 301 of the paper, whereas the 
final dot printed only adds the area 303 between existing dots, as seen in Fig. 3B. 

Letting d represent the area of a single dot, n represent the number of bits of 
colour resolution of a halftoning algorithm, D represent the grid area of a single dot and 
30 f(x) represent the ratio of ink coverage for the domain jc = 0 - 1, such that yj^C^) = 0 and 
f{J) = 7, then: 

Therefore 
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f'(0) « 



I 



D(r-i)' r-i 

D 



= / + p, ; and 



(7) 



for a printer with ideal dot size, 



\r-i) ~'' D(r- 1) ' 



(8) 



10 



Therefore 



fV) 



p-(d.D) I 



D(r-i) ^'r - 1 

D 
d-D 



(9) 



Now, letting 6 = — then f (0) = 1 + S and f (1) = 1 - S. 
If ^is small, as it should be for a good printer, then 



1-5- 



Therefore, 



f'(0) 



f'(l) 



In view of the foregoing, and the fact that a symmetrical section of the function 
15 y Z shows the same characteristic, as shown in Fig. 4 (and indicated as 400), the 



20 



inventor has found that the function = ~ closely fits the characteristics of a single ink 

blend displayed by an ink-jet printer placing ink on paper. With reference to Fig. 4, the 
curve from B to A can be re-scaled to have a range and domain of 0-1 . Therefore, 



- C 



- c 



+ c 



1-c' 



- / 



(10) 



Substituting k for and inverting the function in the range 0 - 1 , 



= 1 



k \_xO 
k \ 1 -k'xd -k) \ 
' I -kY x(l -k) ^- k \ 
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\ k(l 'X)'X(1 -k)^k \ 




(11) 



Cx = (l-f(x))Px ^ f(x)h 



(12) 



The function (11) was substituted in Equation (12) and the result was filled to the 
data for single ink blends using a least squares fit process. The present inventor found 
that the function yielded a good fit, especially to colorimetric data for a cyan blend. 
However, systematic error was found, particularly in the black end of a black ink blend. 

The inventor observed that ink dots have uneven c6lour, possibly caused by 
chromatography. The chromotography was assumed to be responsible for the systematic 
error. To allow for the chromatography effect the inventor based the preferred printer 
characterisation function on a dot 501 with a single halo 502 of different colour, as shown 
in Fig. 5. The inventor also assumed, to simplify the mathematics of the model, that the 
positions of the dot 501 and the single halo 502 were completely independent. The halo 
502 was also allowed to be different in size and intensity for each of the X, Y and Z 
channels. Therefore, the X component of the colour, Cr, of a single ink blend printed on 
a page can be expressed as follows: 



Cx = (I -fix)) (1 'g(x)) Px + f(x) (1 •g(x)) Ix + (1 'f(x)) g(x) Hx + f(x) g(x) HIx (13) 



where: Hx is the X component of the halo colour; and 

HIx is the X component of the overlap of halo and ink. 

The halo size can be assumed to be substantially smaller than the ideal dot size, 
in which case the function g(x) will be of the form: 



where k* is not the same as k. 

When fitting the resulting equation (13) to the blend data, it is not possible to 
determine values for all of the colours Ix and Htx since they are not orthogonal in the 
observation equation set. Therefore, two new constraints are added: 



g(x) = Vx 



(14) 



fix - Ix 
Hix — Ix 

without loss of generality. The equation (13) can be further simplified to 
Cx = (I -f(x)) (1 -k'x) Px + (f(x) + (1 -f(x)) k'x) Ix 



(15) 



(16) 
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The above formula (16) for Cx represents a "Single-Ink Halo Model" and was 
found by the present inventor to show no observable systematic error when blend data 
was fitted to it. 

A least squares approach was taken to deriving the parameters for the simple 
single-ink model (Equation (12)) using the measured observations for X, Y and Z. The 
single-ink model is not linear in k and so it is necessary to derive the parameters 
iteratively as explained below: 

Letting: Px represent the X component of the paper colour; 

Ix represent the X component colour of the ink dot on the 
paper; 

k„^j=kn-^Sk (17) 
PxnH= Pxn+SPx (18) 

Ixn^l = IXn'^Slx (19) 

If koy Pjco and Ixo represent first approximations td the values of Jt, Px and Ix, 
then Equations (17) to (18) represent successive refinements to these approximations, 
provided the three 5 terms can be solved for. 
The observation equations are of the form: 

^ = - x(J.%^J^ ^ '- [x(J-%+^ (20) 

and substituting equations (17) to (19) in Equation (20) results in the following: 

Ignoring second order terms. Equation (21) can be reduced to: 
^""i^' x(J-kJ+k„] x(l-kj + k„ 

■'^^'- xCJ.kJ^kA 
-"^^^[xd-k^ + k^ 

■^Sk[(P^n-Ix.) (,(i%ik^] (22) 
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10 



15 



20 



25 



Equation (22) is now a set of observation equations that can be solved for the 5 
terms using least squares. The /?+/ values are substituted for the n values, and solved 
iteratively until the 6 terms become sufficiently small. 

The same technique of iterative least squares can be employed to solve the Halo 
model (ie. Equation (16)) and the observation equations for the Halo Model are: 



- 1. 



(23) 



Multi-ink characterisation can be modelled by an extension of the above 
Equation (16). As discussed above, the single ink model assumes a mix of two "inks" 
being the ink itself and the ink's halo. For the multi-ink case, the Double Neugebauer set 
of primaries is used. 

By letting S be the set of all inks, then the usual Neugebauer set is the set of 
colours of all combinations of the elements of set S. The set which is similarly defined 
but is based on a set S which contains not only the inks but the halos of these inks as well. 
The basis set S for this set contains double the number of elements. Therefore, the 
number of elements in this set is actually the square of the number in the usual 
Neugebauer set and can thus be referred to as the "Double Neugebauer Set*'. 

The colour of a patch of n inks with levels (a/, ...,«„) for each of the inks is then 
predicted using the "Multi-ink Model" as follows: 



Ylli ^P? Ma,)A-fi(a.,) 
Li=0 



. i=0 



(24) 



where: 



K 



is the Double Neugebauer set of colours 
= {P,,...,P22o }; 
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Cp is the X, Y or Z component of the XYZ colour of the primary p; 

// e p is true if the i'** ink is included in the Neugebauer primary p; 

Hi e p is true if the i*** ink's halo is included in the Neugebauer primary p; and 

gi(ai) = ki'an and 
the expression P? a : b takes the value a if p is true and b if pis false. 
Similarly, the Y and Z components can be computed. 

The multi-ink model (ie. Equation (24)) cannot be solved for k'i and Cpx using 
a minimising method, such as least squares fit, since the least squares matrix is singular. 
The singularity occurs since there is insufficiem information to determine the constants ki 
and /:'/, and the Halo components of the set K. As at Equation (13), the Halo components 
cannot be computed and can be set to the colour without the Halo without loss of 
generality. In addition, some of the multiple Halo components do not contribute 
significantly to the overall result and can therefore be assumed to be either redundant or 
the same as other components. Assuming some of the Halo components to be the same as 
other components significantly reduces the amount of measurement and computation 
necessary to solve the multi-ink model. 

Equation (24) can be solved in two steps, as seen in the flow chart of Fig. 6. The 
process begins at step 601, where the Equation (24) is used to simultaneously solve for 
the single-ink blend parameters k, k\, the paper colour and after inserting the measured 
observation values for X, Y and Z for each of the Neugebauer primaries into 
Equation (16). The solution is non-linear in and so must be carried out iteratively. 
However, it is not possible to solve for and HJx since it is not possible to know all of 
the elements in the Double Neugebauer set. At the next step 603, the remaining multi-ink 
blend parameters are simultaneously directly solved for, using an independent procedure, 
by removing any redundant elements of the Double Neugebauer set. The process 
concludes at step 605, when steps 601 and 603 have been repeated for each of the X, Y 
and Z components. 

As discussed above, the inventor found that using the method of Fig. 6 to solve 
the multi-ink model resulted in the k*i and many of the Double Neugebauer set having 
little physical meaning due to the arbitrary constraints which are imposed at 
Equation (15). In a preferred implementation, step 603 of Fig. 6 requires specific rules in 



527047.doc 



- 14- 



order to decide which of the Neugebauer primaries are redundant and can be eliminated 
from the computation of the multi-ink parameters. Those rules are as follows: 

(i) The values that are already known from the single ink blend computation do 
not need to be calculated again. These known values are the paper colour, single ink 
colours, single halo colours and the combinations of a single ink's colour and 
corresponding halo colour, all of which have been computed in step 601. 

(ii) Where an element includes both colour and halo for any ink, the element is 
assigned to be equal to the element which has the halo removed. 

(iii) Where an element includes two or more halos, the element is assigned to be 
equal to the element which has the halos removed. 

(iv) Where an element includes three or more main ink colours, the element is 
assigned to be equal to the element which has all the halos removed. 

(v) If a combination of inks will flood when printed, there is no requirement to 
know the value of tlie combination accurately. If the combination contains black, it is 
assinged to be equal to black. Otherwise, it is assigned to be equal to the element which 
has all thin inks (ie, C and M' inks) removed. 

The above rules (i) to (v) are applied in sequence to each element of the Double 
Neugebauer set K as the multi-ink parameters are calculated in step 603 for each of the X, 
Y and Z components. 

Three specific sample sets of colour patches, refered to herein as "Measurement 
Sets A, B and C", respectively, are used in the preferred implementation. 

(i) Measurement set A 

Measurement set A comprises ink patches of single ink only. Set A is only used 
in step 601 of Fig. 6 when determining the parameters of single ink blends. For accurate 
printer characterisation, it is necessary to determine the paper colour plus three 
parameters for each ink. If a redundancy of 100% is required to ensure reliable 
determination, at least six sample colour patches for each ink are preferably used. The 
patches should be visually spread between the colour of paper and full ink for the 
particular printer being characterised. 

(ii) Measurement set B 

For accurate characterisation, it is important to have measurements that are 
spread over the whole printer colour gamut to be modelled. However, for most printers, 
there is no real necessity to accurately model the whole gamut since the whole gamut is 
not used in the production of images. 
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Measurement Set B is used to provide information which keeps the multi-ink 
model stable. When used in the least squares process the values measured from set B are 
weighted with a very low weight. 

Measurement Set B comprises all multiple ink Neugebauer primaries which do 
not flood. For a six ink printer there are 64 primaries. Depending on the printer, about a 
third of these will flood and seven are not relevant because they consist of a single or no 
ink. In this circumstance. Measurement Set B will contain 35 patches. For printers with 
fewer than 6 inks, the inventor found that a more stable characterisation was achieved 
when set B was augmented with some mid range samples. 

Set B should hot contain any patches printed with a single ink as these do not 
contribute to the second part (i.e step 603) of the multi-ink characterisation. For a three 
ink printer. Set B is not required since set B will not differ from Set C below, 
(iii) Measurement Set C 

Measurement Set C preferably comprises double the number of samples as the 
number of elements that are required for Measurement Set B. Set C should not contain 
any patches printed with just a single ink as these will not contribute to the second part 
(ie. step 603) of the multi-ink characterisation. Likewise paper colour is not relevant. 

The region of a model gamut most important to accurate printer characterisation 
is that region which results from the "Black Channel Generation (BG), Under Colour 
Removal (UCR) and Ink Split methods" to be used. Measurement Set C comprises of 
colours that are visually evenly spaced throughout this region. 

Table 1 below lists typical sample set sizes and number of parameters that will 
result if the preferred method of characterising a printer is used. Table 1 below lists these 
values for three different ink systems. 

System 1 : CMY inks; 
System 2: CMYK inks; and 
System 3: CVM^'YK inks. 



Table 1 



SYSTEM 


1 


2 


3 


Typical Set A Size 
(Number of Samples) 


7x3 


7x4 


7x6 


Typical Set B Size 


Redundant for 3 inks. 


64 


35 
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(Number of Samples) 




(Includes mid 
range samples) 


(Non-flooding 
Neugebauer 
primaries only) 


Typical Set C Size 
(Number of Samples) 


150 

(Includes mid range 


200 


300 


Step 601 degrees of 
freedom. 


10 


13 


19 


Step 603 degrees of 
freedom. 


13 


34 


125 



The values in Table 1 in relation to the degrees of freedom for steps 601 and 603, 
represent the number of Double Neugebauer primaries required to solve the multi-ink 
model for the specified printer system. As can be observed in Table 1, for the six ink 
printer case, the number of Double Neugebauer primaries required to solve the multi-ink 
model is reduced from 4096 to 125. Therefore, computation is greatly reduced using the 
preferred method of characterising a printer in accordance with the present invention. 

Any printer characterisation model is not useful unless the model can be inverted 
and ink amounts needed to achieve required colours calculated. The Multi-ink model 
(Equation (1)) cannot be analytically inverted. Therefore, a numerical inversion must be 
carried out. Further, any printer characterisation is a function of n variables which 
produces three outputs where n is the number of inks. When /i > J, the function cannot 
be inverted without fiirther constraints. Black Channel Generation (BG), Under Colour 
Removal (UCR) and Ink Split (IS) methods become those constraints. 

Letting C(inks) represent the printer characterisation function, consider the 
function U(c,m,y) which computes the inks based on some arbitrary cyan, magenta and 
yellow values. U(c,m,y) will typically involve BG, UCR and IS. The compound function 
C(U(c,m,y)) can be numerically inverted provided that the function U is chosen 
appropriately. Few generally reliable methods are available for multi-variate root 
determination. In the preferred implementation, multi-variate Newton's method was used 
to solve the equation 

(x^y.z) = C(U(c^m.y)) , 
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for values of (x,y,z) that appear in the printer characterisation gamut. Newton's method is 
iterative and dependent on a first approximation. The inventor found that the choice of 
the first approximation and the choice of the U(c,m,y) function are important in achieving 
a reliable inversion. 

The inventor also observed that on some printers adding yellow ink to black makes 
a lighter colour, presumably because there is a stronger reflective component to yellow. 
Many printers place yellow over black. Therefore, the preferred UCR selected replaced 
100% of yellow but only 80% of Cyan and Magenta. The values of the function U are 
decided for the eight comers of a unit cube. 

U(0,0,0) = (0,0,0,0) 

U(0,0,l) = (0,0,1,0) 

U(0,1,0) = (0,1,0,0) 

U(0,1,1) = (0,1,1,0) 

U(1,0,0) = (1,0,0.0) 

U(1AI) = (1,0,1,0) 

U(1,1,0) = (1,1,0,0) 

U(l,l.l) = (0.2,0.2,0,1) 
Tri-linear interpolation on the resulting ctibe of values is preferably used to compute 
U(c,m,y). Tri-linear interpolation results in a smooth function and Newton's method 
requires a function with continuous first partial derivatives in order to be stable. 

Given the above, a first approximation of (c,m,y) = (0,0,0), that is paper colour, was 
found to be a good approximation for about 80% of in-gamut colours. However, the 
inventor also found that the chosen first approximation was not appropriate for colours 
close to comers of the unit (c,m,y) cube. To improve the reliability, the intersection (in 
XYZ colour space) of the vector through the required colour and printer black with the 
planes (i) to (vi) below is preferably calculated. 



(i) 


C(U(c,m,y)) where (c.m.y) is 


(0,0,0) 


(1.0,0) 


(0.1,0) 


00 


C(U(c,m.y)) where (c,m,y) is 


(0,0,0) 


(0,1.0) 


(0,0.1) 


(iii) 


C(U(c,m,y)) where (c,m,y) is 


(0.0,0) 


(0.0,1) 


(1.0.0) 


(iv) 


C(U(c,m,y)) where (c.m.y) is 


(1,1.0) 


(1.0,0) 


(0,1,0) 


(V) 


C(U(c,m,y)) where (c.m.y) is 


(0,1,1) 


(0,1.0) 


(0,0.1) 


(vi) 


• C(U(c,m,y)} where (c.m.y) is 


(1,0,1) 


(0.0.1) 


(1.0.0) 



The computation is preferably carried out in the order (ie. (i) to (vi)) given 
above and the first computation that determines that the intersection lies within the 
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triangle which defines the plane is preferably used to compute the first approximation. 
The approximation is preferably determined with weighted averaging of the points such 
that the weights are the same as those which result in the intersection point computed. 
The preferred intersection calculation is not possible if the required (X,Y,Z) value is close 
to printer black. Where no appropriate intersection point can be found in this manner, 
white is preferably used. 

If white does not result in a solution, the components of the first approximation 
are preferably halved and Newton's method retried. If still no solution is found (c,m,y) = 
(1,1,1) is used. 

The function C(U(c,m.y)) is not computable outside the domain 0 < c,m,y < 1, as 
any value outside this domain is outside the domain of observations used in the least 
squares process. Therefore, when computing Newton's method, the consecutively refined 
estimates are preferably constrained to the domain 0 < c,m,y < 1 . By using consecutively 
refined estimates as described, the multi-ink model (i.e.: formula (1)) can therefore be 
inverted. 

The above indicates that a method of detemlining a first approximation for the 
iterative computation of amounts of ink which will result in a required colour when 
printed on a printer can be determnined. Where the required colour can be determined by. 
the use of an analytical printer characterisation function, the method operates in a 
Cartesian colour coordinate colour space and commences by determining the colour of 
black on the printer. If the required colour is close to black, the first approximation 
comprises those inks required to produce black. 

If the required colour is not black, then a line is constructed that passses through 
black and the required colour. Then the intersection of the line and each plane defined by 
sets of three colours (cl, c2, c3) is determined, in which each set creating a corresponding 
triangle resulting from the following ink combinations: 



Inkl 


Ink2 


Ink 3 


Paper 


Yellow 


Cyan 


Paper 


Cyan 


Magenta 


Paper 


Magenta 


Yellow 


Yellow + Cyan 


Yellow 


Cyan 


Cyan + Magenta 


Cyan 


Magenta 


Magenta + Yellow 


Magenta 


Yellow. 
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If the intersection of the line and the plane lies within the triangle then it is 
necessary to compute a number of weights (wl , w2, w3) which have to satisfy: 

wlcl + w2c2 + w3c3 = colour at the intersection, 
and then compute the first approximation for the amounts of ink as: 

wl*Inkl + w2*Ink2 + w31nk3, 
otherwise paper colour (no ink) is used at the first approximation. 

In view of the above, the aforementioned preferred method of characterising a 
colour printer provides a relatively easier method of characterising a colour printer using 
a relatively small number of measured samples compared to prior art methods. 

The aforementioned preferred method of characterising a colour printer 
comprises a particular control flow. There are many other variants of the preferred 
method which use different control flows without departing from the spirit or scope of the 
invention. Furthermore one or more of the steps of the preferred method can be 
performed in parallel rather sequentially. 

The method of characterising a colour printer is preferably practiced using a 
conventional general-purpose computer system 700, such as that shown in Fig. 7, wherein 
step 105 of Fig. 1, and the subsequent calculation of the CIE X, Y and Z components, can 
be implemented as software, such as an application program executing within the 
computer system 700. In particular, step 105 is effected by instructions in the software 
that are carried out by the computer. The software can be divided into two separate parts; 
one part for carrying out the method of characterising the printer; arid another part to 
manage the user interface between the latter and the user. The software may be stored in 
a computer readable medium, including the storage devices described below, for example. 
The software is loaded into the computer from the computer readable medium, and then 
executed by the computer. A computer readable medium having such software or 
computer program recorded on it is a computer program product. The use of the 
computer program product in the computer preferably effects an advantageous apparatus 
for determining pixel edge orientation for a pixel-based image in accordance with the 
embodiments of the invention. 

The computer system 700 comprises a computer module 701, input devices such 
as a keyboard 702 and mouse 703, output devices including a printer 715 and a display 
device 714. A Modulator-Demodulator (Modem) transceiver device 716 is used by the 
computer module 701 for communicating to and from a communications network 720, for 
example connectable via a telephone line 721 or other functional medium. The 
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modem 716 can be used to obtain access to the Internet, and other network systems, such 
as a Local Area Network (LAN) or a Wide Area Network (WAN). 

The computer module 701 typically includes at least one processor unit 705, a 
memory unit 706, for example formed from semiconductor random access memory 
5 (RAM) and read only memory (ROM), input/output (I/O) interfaces including a video 
interface 707, and an I/O interface 713 for the keyboard 702 and mouse 703 and 
optionally a joystick (not illustrated), and an interface 708 for the modem 716. A storage 
device 709 is provided and typically includes a hard disk drive 710 and a floppy disk 
drive 711. A magnetic tape drive (not illustrated) may also be used. A CD-ROM 

10 drive 712 is typically provided as a non-volatile source of data. The components 705 
to 713 of the computer module 701, typically communicate via an intercomiected bus 704 
and in a marmer which results in a conventional mode of operation of the computer 
system 700 known to those in the relevant art. Examples of computers on which the 
embodiments can be practised include IBM-PC's and compatibles. Sun Sparcstations or 

15 alike computer systems evolved therefrom. 

Typically, the application program of the preferred embodiment is resident on 
the hard disk drive 710 and read and controlled in its execution by the processor 705. 
Intermediate storage of the program and any data fetched from the network 720 may be 
accomplished using the semiconductor memory 706, possibly in concert with the hard 

20 disk drive 710. In some instances, the application program may be supplied to the user 
encoded on a CD-ROM or floppy disk and read via the corresponding drive 712 or 71 1, 
or alternatively may be read by the user from the network 720 via the modem device 716. 
Still further, the software can also be loaded into the computer system 700 from other 
computer readable medium including magnetic tape, a ROM or integrated circuit, a 

25 magneto-optical disk, a radio or infra-red transmission channel between the computer 
module 701 and another device, a computer readable card such as a PCMCIA card, and 
the Internet and Intranets including email transmissions and infomiation recorded on 
websites and the like. The foregoing is merely exemplary of relevant computer readable 
mediums. Other computer readable mediums may be practiced without departing from 

30 the scope and spirit of the invention. 

The method of characterisation can alternatively be implemented in dedicated 
hardware such as one or more integrated circuits performing the functions or sub 
functions of Fig. 1 . Such dedicated hardware can include graphic processors, digital 
signal processors, or one or more microprocessors and associated memories. 
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The foregoing describes only some embodiments of the present invention, and 
modifications and/or changes can be made thereto without departing from the scope and 
spirit of the invention, the embodiment(s) being illustrative and not restrictive. 

In the context of this specification, the word "comprising" means "including 
principally but not necessarily solely" or "having" or "including" and not "consisting only 
or\ Variations of the word comprising, such as "comprise" and "comprises" have 
corresponding meanings. 
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Tbe claims deflning the invention are as follows: 

1 . A method of predicting colours resulting from using one ink on a printer, said 
method comprising the steps of: 

s producing a test page for said printer, said test page comprising a plurality of 

sample ink patches, each said patch comprising a plurality of ink dots according to a 
predetermined ratio for each patch, said predetermined ratios being based on a number of 
dots printed for a corresponding patch compared to a maximum number of dots; 

measuring each of said sample ink patches to produce a plurality of colour space 
10 coordinates for each said patch; and 

fitting values representing a ratio of dots for each said patch, and said colour 
space coordinates, to a plurality of predetermined functions using a minimisation process 
in order to predict said colours, wherein said predetermined functions are of the form 

P. 0'fx(a)) + hfM 
15 Py(l'fy(a)) + Iyfy(a)'. 

Pz 0'fz(a)) + hf,(ay, 
where a represents a ratio of dots for a particular pzXc\i^fx(ci)^fy(a) mAftfa) are functions 
of a and are derived from a section of the curve y = //c, and wherein Pjr, Py and Pz are 
colour space components representing a paper colour, and /j^, ly and Jz are colour space 
20 components representing the ink colour, wherein said functions fx(o)yfy(Q) and fzia) are of 
the form: 

and wherein a is a variable and is a constant, which can be different in each of fx(cL)^ 
fy(a) and fz(a). 

25 

2. The method according to claim 1, wherein a is equal to zero. 

3. The method according to any one of claims 1 or 2, wherein said functions fx(a), 
fy(a) and fz(a) are derived from an assumption that each ink dot consists of a first colour 

30 being a main ink colour and a second colour being a different colour to said main ink 
colour. 

The method according to claim 3, wherein a is of the form: 
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and wherein k and k' are constants which can be different in each of fx(aj^fy(a) and fz(a). 

5 

5. The method according to any one of claims 1 to 4, wherein said minimization 
process is a least squares process. 

6. An apparatus for predicting colours resulting from using one ink on a printer, 
10 said apparatus comprising: 

means for producing a test page for said printer, said test page comprising a 
plurality of sample ink patches, each said patch comprising a plurality of ink dots 
according to a predetermined ratio for each patch, said predetermined ratios being based 
on a number of dots printed for a corresponding patch compared to a maximum number 
15 of dots; 

means for measuring each of said sample ink patches to produce a plurality of 
colour space coordinates for each said patch; and 

means for fitting values representing a number of dots for each said patch and 
said colour space coordinates to a plurality of predetermined functions using a 
20 minimisation process in order to predict said colours, wherein said predetermined 
functions are of the form 

P. (I 'fx(a)) + //a;; 
Py 0 -fy(a)) + fy(ay. 
Pz a -Ma)) + /z Ma)- 
25 where a represents a ratio of dots for a particular psiichy fx(a)y fy(a) and fz(a) are 

functions of a and are derived from a section of the curve y = J/x^ and wherein P^y Py and 
Pz are colour space components representing a paper colour, and /x, ly and Jz are colour 
space components representing the ink colour wherein said functions fx(a),fy(a) and Ma) 
are of the form: 

and wherein a is a variable and A is a constant, which can be different in each ofMa)y 
Ma) and/;fa;. 
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The method according to claim 6, wherein a is equal to zero. 



8. The apparatus according to claim 6 wherein said functions fx(a)^ fy(a) and fz(a) 
are derived from an assumption that each ink dot consists of a first colour being a main 
ink colour and a second colour being a different colour to said main ink colour. 



10 



The method according to claim 8, wherein a is of the form 



*'"(^ - a(l -% + a) 



and wherein k and k' are constants which can be different in each of fx(ci)yfy(a) and fj(a), 

10. The apparatus according to any one of claims 6 to 9, wherein said minimization 
15 process is a least squares process. 




11. A program for predicting colours resulting from using one ink on a printer, said 

program comprising: 

code for measuring each of a plurality of sample ink patches of a test page for 
20 said printer to produce a plurality of colour space coordinates for each said patch, each 

said patch comprising a plurality of ink dots according to a predetermined ratio for each 

patch, said predetermined ratios being based on a number of dots printed for a 

corresponding patch compared to a maximum number of dots; and 

code for fitting values representing a number of dots for each said patch and said 
25 colour space coordinates to a plurality of predetermined functions using a minimisation 

process in order to predict said colours, wherein said predetermined functions are of the 

form 

Px (1 -Ua)) A h Ua) \ 
Py (1 -fy(a)) + ly fy(ay. 
30 P^ (1 -Ma)) + /, f,(a); 

where a represents a ratio of dots for a particular patch, fj,(aj,fy{a) and f/a) are 
functions of a and are derived from a section of the curve y = 7/r, and wherein P^^ Py and 
Pz are colour space components representing a paper colour, and Ix, ly and Iz are colour 
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space components representing the ink colour, wherein said functions fx(o)^f^o) and fj(Q) 
are of the form: 

a(l X) ^ 

and wherein a is a variable and A: is a constant, which can be different in each of fxio)-, 
5 fy(a) 2x\Af^(a). 



12. The method according to claim 11, wherein a is equal to zero. 

13. The program according to claim 11, wherein said functions fx(a)^fy(a) and f^(a) 
are derived from an assumption that each ink dot consists of a first colour being a main 
ink colour and a second colour being a different colour to said main ink colour. 

14. The method according to claim 13, wherein a is of the form 

and wherein k and are constants which can be different in each of fx(a)yfy(a) and fz(a). 

15. The program according to any one of claims 11 to 14, wherein said minimization 
process is a least squares process. 

16. A method of predicting colours resulting from using one ink on a printer, 
25 substantially as hereinbefore described with reference to Figs. 2 to 6. 

17. An apparatus for predicting colours resulting from using one ink on a printer, 
substantially as hereinbefore described with reference to Figs. 2 to 6. 
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18. A program for predicting colours resulting from using one ink on a printer, 
substantially as hereinbefore described with reference to Figs. 2 to 6. 



DATED this Thirty First Day of October 2002 
Canon Kabushiki Kaisha 
Patent Attorneys for the Applicant 
SPRUSON & FERGUSON 
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